Single exciton spectroscopy of semimagnetic quantum dots 
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A photo-excited 11- VI semiconductor nanocrystal doped with a few Mn spins is considered. The 
effects of spin-exciton interactions and the resulting multi-spin correlations on the photolumines- 
cence are calculated by numerical diagonalization of the Hamiltonian, including exchange interaction 
between electrons, holes and Mn spins, as well as spin-orbit interaction. The results provide a uni- 
fied description of recent experiments of photoluminesnce of dots with one and many Mn atoms as 
well as optically induced ferromagnetism in semimagnetic nanocrystals. 
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Control of the wave function of spins embedded in 
semiconducting materials is in the road-map of the quan- 
tum hardware development"'^. Optical excitation and 
probing of single exciton confined in a nanometric re- 
gion permits to manipulate the quantum state of the 
exciton spin^-'^, or the spin degrees of freedom cou- 
pled to the excitons like nuclei^, donor electrons^S^^ or 
ion,a5,7,8,9,io,ii,i2,i3,i4^ -^Qceni reportsLSia of opti- 
cal detection of the photoluminescence (PL) of a single 
nanocrystal of 11- VI semiconductor doped with a single 
Mn atom showed evidence of a one on one correspon- 
dence between the energy of the emitted photons and 
the quantum state of the Mn spin after photon emission. 
Substitutional Mn in (II,Mn)Vl semiconductors is a neu- 
tral impurity^''' with 5 electrons in the open d shell that 
behave like a spin S — ^. Therefore, the experiment of 
ref. (0) is a proof of principle of the optical manipula- 
tion and detection of the quantum state of a single spin 
5" = 5/2. 

Several other groups have reported the fabrication and 
optical spectroscopy of single nanocrystals of II- VI semi- 
conductors doped with tens of Mn atom3ii*i2iiiiiii&. 
The PL spectrum of one of these semimagnetic nanocrys- 
tal is broader than that of a pure doli^ii^, and it 
shows strong sensitivity to the application of magnetic 
fieldsA^i^, as a result of the exchange interaction of con- 
duction band (CB) electrons and valence band (VB) holes 
with the Mn spins. Here, a quantum theory of the sin- 
gle exciton spectroscopy of a nanocrystal with a few Mn 
spins is presented. The goal is to provide a unified de- 
scription of the PL spectra of nanocrystals with on o^i^i^ 
and many Mn atoms, with emphasis on how to extract in- 
formation about the quantum state of the Mn spins from 
the PL spectra. With that aim, the standard Hamil- 
tonian for semimagnetic nanocrystals is solved exactly 
for Nmu =1,2,3 and 4 Mn ions both when one exci- 
ton is present (excited state manifold, XMS) and absent 
( ground state manifold, GSM). The optical transition 
rates between the GSM and XSM are calculated and PL 
spectra are obtained both in the Faraday and Voigt con- 
figurations for a range of situations. Our results agree 
very well with the PL experiments with dots with one 
Mn£&S and manyiSiH. Previous theory work reporting 



exact diagonalizationsi^ii^ have addressed single Mn dots 
with spherical shape, resulting in PL spectra different 
from those obtained here and reported experimentallji^. 

Hamiltonian. The standard Hamiltoniani^4iLi24i2iS£ 
for diluted magnetic semiconductors describes CB elec- 
trons and VB holes interacting with localized Mn spins 
{Mi) via a local exchange interaction and coupling to 
an external magnetic field, B. For the results presented 
here the short range antiferromagnetic superexchange is 
not effective^" . Both CB electrons and VB holes are con- 
fined in all the 3 spacial directions resulting in a dis- 
crete single particle spectrum. The Hamiltonian reads 
H = Ho + Hi+ Heh, where Ho — -gfisB J2i and 

71,(7.(7' 

^'iiB)dld, + [jeSeiri) + JhSh{ri)) (1) 

In the absence of electron- hole {eh) pairs, the Mi interact 
only with B. The first term in Tii describes the CB 
electrons confined in the nanocrystal orbital levels i'n{r) 
and coupled to the external magnetic field. The Pauli 
matrices are denoted by fa,a' and cjj^ stands for the CB 
electron creation operator. The second term describes 
the VB holes confined in the nanocrystal, including the 
spin orbit interactionSiiS^; the operator dl, creates a VB 
hole in the spin-orbital 0j/(r). The third and fourth term 
describe the exchange between Mn and CB electrons and 
VB holes, respectively. Both the CB-Mn and the VB-Mn 
exchange are local and couple the carrier spin density at 
the Mn location Se,h{ri) to A//. The eh exchange is 
described by Heh — —JehSe ■ Sh- Direct eh Coulomb 
interaction and orbital magnetism are irrelevant in the 
small dots considered hereSS. 

We consider dots whose single particle energy spacing 
is much larger than all the other intra-band energy scales 
of the Hamiltonian. The shape of the dot determines the 
degree of light hole (LH) heavy hole (HH) mixing of the 
VB single particle states which, in turn, determines the 
relative importance of the Ising and spin flip part of the 
hole-Mn and hole-electron exchange interactions. A hard 
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wall cubic potential with dimensions Lx,Ly,Lz conve- 
niently describes the interplay between the shape of the 
dot and the single-hole statesSi, although real dots prob- 
ably have more complicated geometries. We always take 
Lz < Lx, Ly and wc consider light propagation along the 
z axis. The CB electron envelope functions are given by-^ 

V's = ^J^Sin{kxx)Sin{kyy)Sin{kzz) where ki = ^j^, 

i = x,y, z, V = LxLyLz and rii are positive integers. 
The description of the confined VB single particle states 
follows that of reference l|2lh : a zero dimensional ver- 
sion of the 6 band Kohn-Luttinger HamiltonianS^. The 
envelope states for the single particle lowest energy dou- 
blet are V'o(^) = "01,1,1 (^)' both for VB holes and CB 
electrons. The wave functions of the VB lowest energy 
doublet are \+) = ci||, ^) + calf, ^) and |-) ^ K:\+) 
where IC is the time reversal operation and | J, Af) are the 
standard^-^ light hole (LH) and heavy hole (HH) states 
at the Fg point. 

If Lx = Ly LH-HH mixing vanishes (c2 = 0) and 
so they do the matrix elements (±|S'j^'^|=p), resulting in 
a pure Ising (zz) hole-Mn and eh exchange coupling. 
For < \Lx — Ly\ < Lz there is some HH-LH mix- 
ing (0 < |c2| < |ci|) that, on top of the dominant 
Ising coupling, permits some spin flip {xx, yy) hole- 
Mn and eh interaction. In contrast with the VB holes, 
the CB electron-Mn interaction is perfectly isotropio^. 
The results shown here belong to 4 different dots with 
{Lx,Ly,Lz) (7,7,4) nm (dot #1), (15,15,3) nm (dot #2), 
(7,5,2) nm (dot #3) and (7,3.7,2) (dot #4). The strength 
of the exchange coupling between the carriers and the Mn 
spins is je.,h{I) = Je,h\4'o{'i^i)\'^ ■ Foi a given dot the mean 
value of the exchange coupling is {je,h) — -y^ and the 
maximum value is 8{jeM)- For CdTe we have^^ A=+60 
meV nm'^ and Je = — 15meV nm^. We take ge — —1.5 
and = —0.1^. In order to establish the influence of 
the different exchange couplings in the PL spectra of dots 
with 1 Mn, we first take Jeh = (fig. 1) and then we 
study the effect of eh exchange (fig. 2a,b,c) 

The many body states are classified according to the 
number of excitons {Nx — for the GSM and Nx — 1 
for the XSM) and to the number of Mn ions, Nj^jn- The 
number of states in the GSM of a given nanocrystal is 
(6)^Af n . If only the lowest energy doublet of both the CB 
and VB states are kept, there are four possible exciton 
states, labeled by | ± l)x and | ±2)x following their total 
angular momentum, and the number of states in the XSM 
is 4 X (6)^*^". The many-body states are found by stan- 
dard numerical diagonalization of the Schrodinger equa- 
tions Ho|G) = Eg\G) and (Ho+Hi) \X) = Ex\X). The 
GSM states only feature Mn spin coordinates whereas the 
XSM states include Mn spins, as well as 1 CB electron 
and 1 VB hole that occupy linear combinations of the 
single particle states described above. 

In the absence of magnetic field and superexchange 
all the GSM states are degenerate. In contrast, Ising 
hole-Mn coupling of dot #1 lifts splits the 24 levels in 6 
quartets, corresponding to the six possible relative ori- 



entations between the Mn and the hole spin along the z 
direction. Both electron-Mn and electron-hole exchange 
coupling further split the 6 quartets into 12 doublets. 
Depending on the relative size of J^h and Je the low- 
est energy state of the XSM is either | ± l)x "Xi | T |) 
or I ± 2)x iX) I T §)■ The diagonalizations show that the 
lowest energy states of the XSM of dots with Nm =2,3,4 
describe fully polarized Mn spins along Lz- Hence, the 
exact many-spin wave functions of the Mn coupled to 
zero dimensional exciton feature ferromagnetic correla- 
tions, very much like bulk carrier mediated ferromagnetic 
order— and in agreement with experimental observations 
in photoexcited nanocrystalsi^. 

Transitions from the XSM to the GSM are possible 
via spontaneous emission of a photon with the ade- 
quate energy hu and polarization. We consider circu- 
larly polarized photons. The rate of spontaneous emis- 
sion of a photon of energy hiu = Ex — Eg from the 
state \X) to the state \G) is given by the Weisskopf- 
Wigner expression Tax = ■^:^-j\{G\V^\X)\'^ where 

= J2u,n,a{'^H^ ± «2/)l"'0-)44,<T + h.c. is the sec- 
ond quantization representation of the interband electric 
dipole operator that yields the standard optical selec- 
tion rules. Standard optical selection rules forbid photon 
emission from | ± 2)x states. Since the electric dipole 
operator does not affect the Mn d electrons, the Mn spin 
part of the collective wave function does not change dur- 
ing the recombination process. As a result, most of the 
GSM to XSM transitions that are allowed by the stan- 
dard optical selection rule but are forbidden due to the 
orthogonality of the spin part of GSM and the XSM 
states. This optical spin blockade, reminiscent of the spin 
blockade well characterized in single electron transport-^, 
implies a huge limitation to the otherwise large set of 
XSM to GSM transitions. The PL spectrum is calcu- 
lated by statistical weighting of the initial states \X): 

Ipl{,oj) = PWgxS [hu - {Ex - Eg)) 

X,G 

where P{X) is the occupation probability of the XSM 
state \X). In general the calculation of P{X) would in- 
volve the determination of the non-equilibrium steady 
state density matri5iiii&. The experimental results of 
reference support approximating P{X) as a thermal 
distribution function with an effective temperature. 

In the following we show PL spectra calculated follow- 
ing the above recipe. The spectral lines are broadened 
using gaussians with a FWHM of either 7 =50/ieV or 
7 =400/ieV, adapted to the experimental photodetector 
resolution of references Qj and fl3l respectively. Sin- 
gle spin optical spectroscopy calculations are shown in 
fig. 1. The zero field PL for nanocrystal #1 with a sin- 
gle Mn impurity and 2 effective temperatures ,0.5 and 2 
meV (6 and 24 Kelvin) is shown in fig. la. The model ac- 
counts quantitatively for the experimental observations^: 
only 6 lines are clearly seen, in spite of the fact that the 
XSM has 12 non degenerate lines, that occupy a spec- 




FIG. 1: (Color online), (a) PL for dot #1 with 1 Mn 
atom and two values of T. (b) Magneto PL in the Faraday 
((0,0, 5r)) and Voigt ((5T,0,0)) configurations. (c,d) Fox 
for the Faraday (c) and Vo.igt (d) cases. 



FIG. 2: (Color online). Upper panels: Influence of the LL- 
HH mixing and the eh exchange on the PL of single Mn PL 
for dots #3 (2a) and #4 (2b and 2c). 2d-2g Zero field PL 
spectrum for dot #2 with 3 Mn impurities as a function of 
the relative Mn-carrier couplings. 



tral range of 1 meV. The number of the lines differ from 
those of previous theory work (liStiS) due to the different 
LH-HH mixing of their model. The 6 hnes undergo a po- 
larization dependent splitting (fig. lb, lower spectrum) 
upon application of a magnetic field in the Faraday con- 
figuration, B = {0,0, Bz), in good agreement with the 
experiments ''■ ^'^ . 

The detection of a photon with a given polarization 
and energy in one of these 6 peaks yields information 
about the final state of the GSM in the optical transition. 
That information is given, for the Faraday spectrum of 
fig. lb, in the fig. Ic. There the size of the symbols is 
proportional to Tqx for (t_ (open squares circles) and 
CT+ (solid diamonds). It is seen how each member of the 
XSM is coupled via a given circularly polarized photon 
with, at most, one member of the GSM. This fact permits 
to map the photon state (energy and polarization) to the 
final Mn state and is the second important result of this 
paper. In particular, going from low to high energy the 
six peaks from the PL in the Faraday case correspond to 
final Mn states going from Sz = +^ to Sz = The 
zero field case has the same one-on-one correspondence. 

The calculated PL in the Voigt configuration, B = 
{Bx,0,0), is remarkably different. The shape of the PL 
is modified drastically, and the spectrum is much wider. 
The GSM states are the eigenstates of the Sx spin oper- 
ator, whereas the XSM states have strong overlap with 
the eigenstates of the Sz operator, because the spin of 
the hole is pinned in that direction and the hole-Mn Ising 
coupling is dominant. The matrix {G\V\X) is now pro- 
portional to {Sx\Sz) making spin blockade much less ef- 
ficient, as shown in fig. Id. This is seen in figure Id and 
accounts for the PL shape reported inlTol 

We now consider the effect of the eh exchange interac- 
tion, absent in the results of figure 1. The Ising part 
of the eh interaction splits bright and dark excitons, 
whereas the spin-flip part of the eh interaction, which 
is only possible if there is LH-HH mixing, mixes -fl and 



— 1 excitons resulting into linearly polarized spectra^. In 
the presence of Mn, this mixing competes with the Ising 
hole-Mn coupling. The combined action of electron-Mn 
exchange and the transverse hole-Mn and eh interactions 
mix somewhat bright and dark excitons. In figure 2a-2c 
we show the PL spectra corresponding to dot ^3 (2a) 
and dot #4 (2b, 2c), all of them with a single Mn atom 
and Jeh = 1 meV. Dot ^^3 features some LH-HH mixing. 
Ising eh interaction red shifts the dark excitons, result- 
ing in the appearance of 3 low energy peaks (compared 
with fig. 1). The middle PL corresponds to a Mn weakly 
coupled {jh = O.WmeV, je — 0.25 jh) to the exciton and 
a strong LH-HH mixing. The PL is linearly polarized 
and the peak structure of figs, la is not resolved. The 
lower panel corresponds to the same dot with a strongly 
coupled Mn {jh = 0.31 meV). There the 6 peak structure 
is recovered, but the PL is linearly polarized, mostly in 
the central peaks for which the hole-Mn Ising coupling is 
smaller. These whole picture is consistent and qualita- 
tively identical to the experiments recently published^. 

We finally address the central theme of this paper: how 
the PL spectra of the dot evolve as the number of Mn 
atoms increases. In this regard, we start considering the a 
dot with 3 Mn atoms, two of which are identically coupled 
(j/i(l) = j/i(2) — 0.38 meV) and the third changes, so 
that the relative coupling to the fermions defined by r = 

is varied from from 1 to 0.1 (figs 2d to 2g). The 
symmetric case r = 1 shows a clean PL spectrum with 
16 lines corresponding to quantum states with the 25* -I- 1 
possible orientations of the collective Sz of the 3 Mn 
atoms that behave like a single magnet with {S = 15/2). 
These 16 peaks are still seen for r = 0.85. At r = 0.5 
the PL is a rather featurless wide peak. When r = 0.1 
the PL spectrum features 11 peaks, that are the analog 
of figure 2d for 2 Mn atoms behaving collectively with 
S = 10/2. It is apparent that the weakly coupled Mn 
spin does not affect the PL. 
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The results of figure 2d-2g suggest that a few Mn 
atoms, strongly coupled to the exciton could be respon- 
sible of the main features of the PL signal observed in 
earlier experiments. In ref. (T^) the PL spectrum of sin- 
gle semimagnetic nanocrystals show a zero field cr+ line 
with a full width at half maximum (FWHM) of ~ 5 meV 
that red shifts and narrows upon application of a mag- 
netic field in the Faraday configuration^^-^^'^'^. This has 
motivated statistical interpretation of the PL spectrum 
of single semimagnetic nanocrystalsiSii^, linking the PL 
lincwidth to the magnetic statistical fluctuations. In such 
an approach the PL FWHM has been shown to scale 

with ^ ^f^. On the other side, it is apparent that the 
PL spectrum of a single Mn impurity^ lies in a window 
of order A/ = | {\jh{l)\ + \je{l)\), which scales like V~^. 
The PL linewidth of fig. 2 also scales with A/ oc . 

We have explored whether the quantum approach with 
Nmu =3 and 4 can model PL spectra like those of refer- 
ences ("11 1121131) . The justification lies in the distribution 
of couplings jh(I), shown in fig. 3a for Cdo.99Mno.01Te 
(dot #1) obtained by random generation of 100 realiza- 
tions of the Mn disorder. Such a dot has Nmu =28 and 
{jh) = 0.3 meV. However, 65 % of the Mn have jt < (jh) 
and, in average, the 4 Mn with largest coupling have 
jh > 3(j/i). We calculate the PL taking the 4 Mn atoms 
with the largest overlap to the carrier wave function, for a 
randomly selected Mn disorder realization. Their jh are 
1.69, 1.66, 1.1, and 0.45 meV. In figure 3b we show the 
corresponding (T+ PL spectrum for values of the magnetic 
field going from zero to 12 T in the Faraday configura- 
tion. We take 7 — OAmeV and T = 3meV. As in the 
experiments, the cr+ PL red shifts as goes up. The 
featureless zero field PL spectrum (with a FWHM of 5 
meV) develops, as Bz increases narrow structures, as re- 
ported experimentallyi^. The narrowing is also seen in 
fig. 3c, that shows the PL of dot #1 with 3 Mn atoms, 
with values of jh closer to the average (0.39,0.40 and 0.44 
meV), with T = 2 meV and 7 = 0.4 meV. In fig. 3d we 
show how the corresponding FWHM decreases both upon 
cooling and upon increasing B^ ■ The former is due to the 
reduction of the spectral range of the occupied states in 
the XSM, resulting as well in a red shift thas has been 
observed experimentallyi^. The magnetic field produces 
a similar effect. 

In conclusion, the problem of a few Mn atoms inter- 
acting with a single exciton in a quantum dot is solved 



via numerical diagonalization of the Hamiltonian includ- 
ing electron-Mn, hole-Mn, electron-hole exchange as well 
as spin-orbit interactions and the correponding PL spec- 
tra are calculated. The model accounts for a series of 
recent reports with the PL of CdTe dots doped with a 
single Mn atom and permits to link the quantum state of 
one Mn spin to the detection of a photon with a given en- 
ergy and polarization. The crucial ingredient for this one 
on one correspondence is the anisotropy of the hole-Mn 
interaction, due to strong spin-orbit interaction for the 
holes. The model reproduce qualitatively the PL spectra 
of single dots with presumably tens of Mn atomsiitiSii^ 
with just 3 and 4 Mn spins strongly coupled to the ex- 
citon and with different coupling strengths. This might 
indicate that the width of the PL spectra coming from 
a single semimagnetic dot might arise from a few Mn 
atoms strongly coupled to the exciton rather than from 
statistical fluctuations of many Mn atoms. 
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FIG. 3: (Color online), (a) Normalized histogram of Mn-hole 
coupling. The dashed line is the integrated histogram, (b) PL 
spectrum for 4 Mn impurities (see text) for several values of 
B;,. (c) Idem for 3 Mn impurities . (d) FWHM for as a 
function of and r=3,2,l and 0.5 meV (top to bottom). 



^ B. E. Kane, Nature (London), 393, 133 (1998). D. Loss et 

al, Phys. Rev. A 57, 120 (1998). 
^ N. H. Bonadeo et ai, Science 282,1473 (1998) 
^ Xiaoqin Li et al. Science 301, 809 (2003). 
* F. Jelezko et al, Phys. Rev. Lett. 93, 130501 (2004) 
^ J. Bao, et al, Nature Mater. 2, 175 (2003). 
^ P. Chen et ai, Phys. Rev. B69 , 075230 (2004) 

L. Besombes et al. Phys. Rev. Lett. 93, 207403 (2004) 



* L. Besombes et al. Phys. Rev. B. 71, 161307 (2005) 
^ Y. Leger et al, Phys. Rev. Lett. 95,047403 (2005) 
^° Y. Leger et ai, submitted to PRB 

A. A. Maksimov et al. Phys. Rev. B62, R7767 (2000) 
G. Bacher et ai, Phys. Rev. Lett. 89 , 127201 (2002) 
P. S. Dorozhkin et al., Phys. Rev. B 68, 195313 (2003) 
S. Mackowski et al. Appl. Phys. Lett. 84, 3337 (2004). 
J. K. Furdyna, J. Appl. Phys 64 R29 (1988). 



5 



S. C. Erwin et al, Nature 436, 91 (2005) 

A. Efros, E. Rashba, M. Rosen, Phys. Rev. Lett. 87, 

206601 (2001) 

A. O. Govorov, Phys. Rev. B 70, 035321 (2004) A. O. 
Govorov and A. V. Kalameitsev Phys. Rev. B 71, 035338 
(2005) 

A. K. Bhattacharjee and J. Perez-Conde Phys. Rev. 

B68,045303 (2003) 

J. Fernandez-Rossier and L. Brey, Phys. Rev. Lett. 93 



117201 (2004) 

F. V. Kyrychonko and J. Kossut, Phys. Rev. B 70, 205317 
(2004) 

22 J.M. Luttinger, W. Kohn, Phys. Rev. 97 869 (1955) 

23 H. Ohno, Science 281, 951 (1998). T. Dietl et al, Science 
287 1019 (2000). 

24 J. J. Palacios et. al., Phys. Rev. B50 5760 (1994) 



